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Abstract 

We  consider  the  asymptotic  solution  of  boundary  value  problems  for 
the  vector  system 


x = A ( t , £ ) x + B (t , c) y + C(t,e) 
zy  = E(t,e)x  + F(t,E)y  + G(t,e) 


as  e ->  0 under  the  assumption  that  the  matrix  F(t,0)  is  singular.  A 
full  set  of  asymptotic  solutions  is  constructed  assuming  that  F(t,0)  can 
be  block-diagonalized,  the  reduced  problem  is  consistent,  and  a new  stabil- 
ity condition  holds.  Boundary  value  problems  are  then  solvable  if  an 
appropriate  "boundary"  matrix  is  nonsingular  for  e / 0.  Such  problems 
arise  in  optimal  control  theory,  among  other  applications. 
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A Singular  Singularly-Perturbed  Linear  Boundary  Value  Problem 


1.  Introduction. 

Let  us  consider  a linear  system  of  the  form 

x = A(t,e)x  + B(t,e)y  + C(t,s) 

ey  = E(t,e)x  + F(t,e)y  + G(t,e) 

for  vectors  x and  y of  dimensions  n and  m,  respectively,  for  a small 
positive  parameter  e,  and  for  a finite  t interval,  say  0 < t < 1.  It 
is  natural  to  consider  (1)  subject  to  a list  of  n + m linearly  independent 
boundary  conditions  of  the  form 

1 

(2)  £ (R  (e)x(j)  + S (e)y(j))  = c(e) 

j=0  J J 

and  study  the  asymptotic  solution  of  (1)  - (2)  as  e •+  0. 

We  recall  that  rather  classical  methods  can  be  used  to  asymptotically 
solve  the  "regular"  singularly  perturbed  problem  (1)  - (2)  when  F(t,0) 
satisfies  an  exponential  dichotomy,  i.e.  its  eigenvalues  have  either  a pos- 
itive or  a negative  real  part  throughout  0 £ t < 1 (cf.,  e.g.,  O'Malley 

(1969a),  Harris  (1973),  or  Ferguson  (1975)).  When  F(t,0)  is  everywhere 
stable,  for  example,  they  show  that  the  initial  value  problem  for  (1)  has 
a unique  solution  which  converges  as  e -►  0 for  t > 0 to  the  solution 
of  the  reduced  system 


□ □ 
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(3) 


/ . 

XQ  = A(t,0)XQ  + B(t,0)Y0  + C(t,0) 


0 = E(t,0)XQ  + F(t,0)\’0  + G(t,0) 


subject  to  the  initial  condition  Xq(0)  = x(0).  Kor  the  analogous  terminal 
value  problem,  however,  the  solution  would  then  be  exponentially  large  as 
e -*  0 for  t < 1.  More  generally,  boundary  layers  (regions  of  nonuniform 
convergence)  of  thickness  0(e)  must  be  expected  at  each  endpoint  for  regular 
problems  when  the  limiting  solution  within  (0,1)  is  bounded,  and  this 
limiting  solution  must  satisfy  the  reduced  system  (3)  and  n boundary  condi- 
tions determined  by  an  appropriate  combination  of  the  original  conditions  (2) 
evaluated  at  e = 0.  Because  such  limiting  solutions  involve  only  n bound- 
ary conditions,  m linearly  independent  solutions  of  the  homogeneous  form  of 
(1)  are  of  boundary  layer  type,  i.e.  they  are  asymptotically  negligible  away 
from  the  endpoints.  Considerable  progress,  then,  has  been  made  in  determin- 
ing which  regular  singularly-perturbed  linear  boundary  value  problems  have 
limiting  solutions  as  e -*■  0,  what  boundary  value  problems  these  limits 
satisfy  within  (0,1),  and  the  nature  of  the  endpoint  boundary  layers.  The 
results  are,  however,  more  complicated  than  for  scalar  problems  (cf.  O'Malley 
(1969b)  and  O'Malley  and  Keller  (1968)).  In  addition  to  its  direct  utility, 
that  information  is  useful  in  analyzing  nonlinear  problems  (cf.  Hoppensteadt 
(1971))  and  in  designing  numerical  algorithms  (cf.  Flaherty  and  O'Malley 
(1977)). 

Here,  we  shall  consider  "singular"  problems  where  F(t,0)  is  singular 
and  of  constant  rank  throughout  0 j t < 1.  Their  analysis  and  the  behavior 
of  their  solutions  are  considerably  more  complicated  than  for  regular  prob- 
lems. Specifically,  we  shall  find  that  the  asymptotic  analysis  of  singular 


i 
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problems  involves  a consistency  condition  which  did  not  occur  for  regular 
problems,  a new  stability  requirement,  and  the  occurrence  of  other 
(thicker)  boundary  layer  regions  of  nonuniform  convergence.  Those  singular 
problems  are  less  complicated,  however,  than  turning  point  problems  where 
F(t,0)  is  singular  at  isolated  points  (cf.  Levinson  (1951),  Wasow  (1965), 
and  Olver  (1977)).  Fundamental  matrices  for  homogeneous  systems  (1)  without 
turning  points  can  be  constructed  as  in  Turrittin  (1952),  and  they  could,  in 
theory,  be  used  to  asymptotically  solve  nonhomogeneous  problems  via  variation 
of  parameters. 

Our  interest  in  such  problems  arose  in  analyzing  nearly  singular  optimal 
control  problems  (cf.  O'Malley  and  Jameson  (.1975,  1977))  and  in  devising 
methods  for  the  numerical  integration  of  stiff  differential  equations  (cf. 
Flaherty  and  O'Malley  (1977)).  The  technique  we  use  generalizes  that  devel- 
oped for  singular  arc  computations  (cf.  Goh  (1966)  and  Robbins  (1967)). 

Closely  related  methods  are  given  in  O'Malley  and  Flaherty  (1977)  and  O'Malley 
(1978)  for  certain  nonlinear  problems. 

The  simplest  control  example  is  given  by 

( 

x = -y,  x ( 1 ) = 0 

< 

cy  = -x,  y (0)  = -1. 

Here,  -x/c  represents  an  optimal  control  and  y,  the  corresponding  state 
of  a nearly  singular  control  problem  (cf.  O'Malley  and  Jameson  (1975)).  The 


x(t)  = -/e  (e 


Z//1  - e‘(1-t)//?)/(l 


+ e 


-2/ 


/c. 


solution  is 


and  y = -x,  so  Che  limiting  solution 


(x(t)  ,y  (t) ) * -(*^,l)e  t/,£ 


is  asymptotically  trivial  for  t > 0 and  features  a 0(/e)  boundary  layer 
at  t = 0.  Note,  in  particular,  that  the  corresponding  control  acts  like 
an  initial  delta-function  impulse. 


2 . A Transformed  Problem  and  the  Corresponding  Reduced  System. 

Under  rather  mild  assumptions,  the  singular  matrix  F(t,0)  can  be 
block-diagonalized  (cf.  Sibuya  (1958,  1966)  and  Chapter  VII  of  Wasow  (1965)) 
We  shall  simply  assume 


(Hi)  that  there  exists  a_  smooth  nonsingular  matrix  P(t)  such  that 


P (t)F(t,0)P(t)  = 0 F5(t,0)  0 


F9(t,0) 


where  - F^  and  F^  are  stable  matrices  (i.e.,  their  eigenvalues  have 
negative  real  parts)  throughout  0 <_  t <_  1 oj[  dimensions  m 2 * m^  and 
m3  * m3’  resPectively»  with  m = + m^  + m^,  m^  > 0. 


Hypotheses  guaranteeing  the  existence  of  P are  given  in  Wasow  (1964)  and 
elsewhere.  (We  note  that  an  analogous  trichotomy  was  used  by  Hoppensteadt 
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and  Miranker  (1976),  except  that  they  allowed  F(t,0)  to  have  purely 
imaginary,  but  nonzero,  eigenvalues.)  In  analogy  with  the  parallel  situa- 
tion in  singular  optimal  control  (cf.  Jacobson  (1971)  or  Anderson  (1973)), 
we  might  call  problems  where  m = m^  totally  singular  and  those  where 
m > m^  partially  singular.  (If  either  m2  = 0 or  m^  = 0,  it  may  be 
more  convenient  to  use  a singular  value  decomposition  F = UDV  where  U and 
V are  orthogonal  and  D is  diagonal  with  the  eigenvalues  of  /f' F.  One 
then  uses  the  transformed  variable  z = U'y  (cf.  O'Malley  (1978)). 

In  general,  we  introduce 


(5) 


,-l 


P > = (yj  y'2  y’)' 


in  (1)  and  obtain  the  equivalent  system 


/ x = A(t,e)x  + B1(t,e)y1  + B2(t,e)y2  + B3(t,e)y3  + C(t,e) 


(6)  \ 


eyi  = Ei(t’c)x  + eF1(t»E)y1  + eF2(t,E)y2  + eF3(t,e)y3  + G1(t»E) 


ey2  = E2  ’ G^x  + EE4(t»E)y1  + E5(t,c)y2  + cF6(t,e)y3  + G2(t,e) 


\ ey3  = E3(t,e)x  + eF7(t,e)y1  + £Fg(t,e)y2  + F9(t,t)y3  + G3(t,e) 


compatible  with 

(A), 

-1 

eF1 

eF2 

eF 

P (FP  - eP)  = 

eF4 

F5 

eF 

1 

eF7 

CF8 

F, 

1 


J d 
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BP  = [B1  B2  B3J, 


E = [E*  E’  E^]', 


and 


P 1G 


[G1  G2  C3]'- 


Experience  with  singular  perturbation  problems  leads  us  to  expect  that  the 
limiting  solution  to  (6)  within  (0,1)  will  satisfy  the  reduced  system  ob- 
tained by  setting  e = 0 in  (6),  i.e. 


( X0  A0X0  + B10Y10  + B20Y20  + B30Y30  + °0 


(7)  < 


° F'10X0  + Gio 


° = E20X0  + F50Y20  + G20 


° E30X0  + F90Y30  + G30 


where,  e.g.,  = A(t,0).  For  any  solution  of  (7), 


(8) 


Y20  ~ F50(E20X0  + G20)’ 


V Y30  F90(E30X0  + G30) 


and  there  remains  the  linear  equations 


(9) 


E10x0  -G10 


tH 

and  the  n order  system  of  differential  equations 


xo  = Hoxo  + B10Y10  + Jo 


(10) 
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to  determine  the  m^-vector  Y^q  and  the  n-vector  X^.  Here 

H0  = A0  " B20F50E20  " B30F90E30  and  J0  = C0  ~ B20F50G20  ~ B30F90G30' 
In  order  to  solve  (7),  then,  we'll  assume 


(H2)  G is  in  the  range  of  for  0 < t < 1 


(H3)  _E10B10  — sta'3^-e  throughout  0 _<  t _<  1. 


Differentiating  (9)  we  have 


E10X10  E10X10  " G10’ 


so  (10)  and  (H3)  allow  us  to  solve  for  as 


(H)  Y1()  - (~E10B10)  [(EioH0  + E10^X10  + ^E10J0  + G10^  * 


and  there  remains  a nonhomogeneous  system  for  X^.  To  determine  Xq,  it 
is  convenient  to  introduce  the  projection 


E ■ •„  - «ho 


where 


Q ■ ,10<E10B10>  • 


We  note  that  EinE  = 0 while  Eir.Q  = I and  EBir.  = EQ  = 0.  Moreover, 
10  10  10  9 


W 
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and  have  rank  since  ^as  rank  m^.  Thus,  E 

has  rank  n - _>  0.  Indeed,  Q is  nearly  a generalized  inverse  of  E^q 

(cf.  Campbell,  Meyer,  and  Rose  (1976)).  Using  (9),  (12)  implies  that 


xo  Exo  ~ QG10 


and  (10),  (11),  and  (14)  imply  the  linear  system 


(EXQ)  = K(EXQ)  + L 


for  EXQ  where  K = EHq  - QE-^q  and  L = -KQG^q  + Ej^  + QG^q.  Under  hypo- 
theses (HI)  - (H3) , then,  the  solution  of  the  reduced  system  (7)  will 

be  completely  and  uniquely  determined  up  to  later  specification  of  a 
boundary  value  for  EXg.  It  is  perhaps  most  natural  to  use  the  condition 


E(j)XQ(j)  = E(j)x(j),  j = 0 or  1, 


presuming  a boundary  value  (16)  is  supplied  by  (2).  Other  possibilities 
should  also  be  considered,  however. 

Note  that  our  manipulations  allowed  us  to  determine  E^qXq  from  the 
linear  equation  (9)  and  the  remaining  "component"  EX^  of  Xq  from  an 
end  value  problem  like  (15)  - (16).  The  alternative  problem  character  of 
the  solution  for  Xq  (cf.  Cesari  (1975))  makes  it  quite  different  from  the 
more  straightforward  solution  of  regular  problems. 

If  (H2)  fails,  the  reduced  system  (7)  is  inconsistent,  but  irrelevant 
(cf.  O'Malley  (1978)).  A simplified  example  is  provided  by  ey  = 1, 
y(0)  = 0 where  the  limiting  solution  for  t > 0 is  unbounded  like  t/c. 


For  regular  problems,  the  reduced  problem  (3)  is  necessarily  consistent. 
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but  inconsistency  of  (3)  would  occur  here  if  (H2)  failed.  Although  we  have 
only  used  the  nonsingularity  of  E10B10  t0  deEine  t,ie  reduced  problem,  the 

stability  assumption  (H2)  is  generally  needed  in  order  for  a limiting  solu- 

• • 

tion  to  exist.  A simple  example  is  provided  by  x = -y,  ey  = x,  x(l)  = 0, 
y(0)  = 1 which  has  the  solution  (x,y)  = (-i/e  sin  t/fc,  cos  t/v^)  for 
which  there  is  no  limit  as  e ->  0.  We  note  that  by  changing  the  sign  so 


that  ey  = -x,  E]qB10  < and  we  have  a limiting  solution.  No  such  stab- 


ility assumption  was  required  for  regular  problems  where  = 0.  If 


E = 0,  further  differentiation  of  (9)  might  allow  one  to  determine  Y , 


'10  10  " ‘10’ 

just  as  singular  arcs  of  higher  order  are  obtained  in  control  (cf.,  e.g., 

Robbins  (1967)).  The  structure  of  the  asymptotic  solutions  will  then  differ 

considerably  from  when  ( H 3 ) holds.  An  example,  arising  in  optimal  control,  is 


-yl  + x2 ’ 

x1(l)  = 0 

“Xl» 

x2(l)  = 0 

ey2> 

yx(0)  = 1 

-x2  - cyl5 

y2(0)  = 0 

(cf.  O'Malley  and  Jameson  (1977)).  Here,  the  asymptotic  solution  is  given 

by  x,  = -X£  and  y ^ = y^  where  X2  = 2 /e  Im  [ce  *?F 

„ „ r -wt/  7e1  , 

y^  = 2 Re  [ce  ] for  c = (1  - 1 


and 


- i \[(  1 + i/c)/(l  - i/e))  1 and 

7T 1 / A J J — — 

w = e VI  + i/e.  Thus,  the  boundary  layer  thickness  is  0(  /c).  Finally, 

if  is  singular,  but  nonzero  (as  when  n < in^) , progress  could  be 

generally  made  through  preliminary  algebraic  manipulations  (cf.  Anderson 

(1973)  for  an  analogous  control  problem). 
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3.  Construction  of  Asymptotic  Solutions. 


A linear  nonhomogeneous  boundary  value  problem  can  be  solved  by  varia- 
tion of  parameters  once  a complete  set  of  linearly  independent  solutions  of 
the  corresponding  homogeneous  system  is  known.  For  the  asymptotic  solution 
of  (1),  we'd  need  n + m linearly  independent  asymptotic  solutions.  Alter- 
natively, one  could  seek  an  outer  solution  of  (1)  (i.e.,  a regular  perturba- 
tion of  an  already  obtained  solution  of  the  reduced  system  (7))  and  modify 
it  by  adding  the  appropriate  boundary  layer  corrections  (cf.  O'Malley  (1969b) 
which  solves  corresponding  scalar  problems).  Since  we  can  supply  n - 
boundary  conditions  (like  (16))  for  the  reduced  system,  we  can  expect  an 
outer  solution  under  hypotheses  (Hi)  - (H3)  to  be  of  the  form 


(X ( t , /e) .Y^ft,  Sz) ,Y2(t, /e) ,Y^ (t,  il)) 


* } (Xj(t),Ylj(t),Y2j(t),Y3j(t))eJ 


being  an  asymptotic  solution  within  (0,1)  which  converges  to  a solution 
(Xq.Y^q, y2q,Y3q)  of  the  reduced  system  (7)  as  e * 0.  It  would  be  completely 
determined  termwise  by  the  boundary  values 


E(j)X(j,/e),  j = 0 or  1 


since  higher  order  terms  will  satisfy  a system  of  the  form  (7)  with  succes- 
sively known  nonhomogeneous  terms. 


Since  E has  rank  n - m^ , the  outer  solution  is  parameterized  by 
n - m^  vector  functions  of  /e,  and  there  is  need  for  m + n - (n  - m^) 
m + linearly  independent  boundary  layer  solutions  of  the  homogeneous 
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version  of  (6)  which  are  asymptotically  negligible  within  (0,1).  For  the 
regular  problem  with  m2  + m^  = m,  there  would  be  m^  boundary  layer  solu- 
tions which  are  decaying  functions  of  the  stretched  variable  k = t/e  and 
m2  which  are  decaying  functions  of  p = (1  - t)/e  (cf.,  say,  Harris  (1973)), 
For  our  singular  problems,  however,  we  shall  also  find  m^  boundary  layer 
solutions  depending  on  each  of  the  stretched  variables  t = t//c  and 
a = (1  - t)//e.  (These  thicker  (/e  >>  e)  boundary  layers  (upon  matching) 
now  require  our  asymptotic  expansions  to  all  be  power  series  in 


y = ve, 


rather  than  e.  In  order  to  generate  these  asymptotic  solutions,  we'll  as- 


sume the  coefficients  in  (1)  to  be  infinitely  differentiable,  though  finite 
approximations  could  be  obtained  under  less  smoothness. 

We  shall  construct  asymptotic  solutions  to  (6)  of  the  form 


, . . . a,  , 6+1  , , , y+2  , . , 6+2 . . . 

x(t,e)  = X(t,p)  + y m ( t , p ) + y n(a,p)  + p r(p,y)  + p j(x,p) 

yx(t,e)  = Y1(t,p)  + pa  1p1(t,p)  + pPq1(a,p)  + pY+2s1(p,p)  + p'5+2J’1(k,p) 

y2(t,c)  = Y2(t,y)  + pap2(x,p)  + pf3+1q2(a,p)  + pYs2(p,p)  + p6+2£.2(k,p) 

y3(t,e)  = Y3(t,y)  + yap3(T,y)  + pP+1q3(a,p)  + pY+2s3(p,p)  + p6£3(tc,p) 


where  the  functions  of  the  stretched  variables 


(21)  x = t/p,  o = (1  - t)/p,  p = (1  - t)/p2,  and  k = t/p2 


tend  to  zero  as  the  appropriate  variable  tends  to  infinity.  The  outer  solu- 
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tion  (17),  then,  provides  the  asymptotic  solution  within  (0,1).  The 

scalings  pQ,  pB,  p^ , and  p^  for  the  boundary  layer  corrections  remain 

free  to  meet  various  boundary  conditions  (2),  while  the  remaining  p \ 

2 

p,  and  p factors  are  simply  used  to  prevent  calculation  of  trivial 
coefficients . 

Since  the  outer  solution  (17)  satisfies  the  nonhomogeneous  system  (6), 
the  boundary  layer  corrections  will  satisfy  the  corresponding  homogeneous 
system.  The  boundary  layer  solutions  which  depend  on  k or  p are  com- 
pletely analogous  to  those  obtained  for  the  regular  problem  when  = 0. 

Let  us  first  seek  m^  boundary  layer  solutions  of  the  homogeneous 
system  (6)  of  the  form 


1 lk 

(m>  “ PI’P2'P3)  ^ 1 ^ plk(x)  ,p2k(T)  ,p3k(T))lJ 

lc  0 


(cf.  (20)).  Thus,  we'll  have 


- ~ = Am  + - B p + B0p„  + B„p_ 
p di  p 11  2‘  2 3r3 


dT  = E!m  + ,,F1P1  + p2p2P2  + p2p3P3 


p ~ = E2m  + pFaPi  + F5p2  + p F6p3 


p — = E3m  + pF7Pi  + p F8p2  + F9p3. 


When  p = 0,  this  reduces  to  the  limiting  problem 


dx  B10(0)plC 
dp10 

-dT  = Eio(0)mo 
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(24) 


( 0 ‘ E20(0>”0  + F50(0)>,20 


0 - E30(0)n,0  + F90(0)p30 


SO 


(25) 


P20(T) 


p30^T^ 


■F-J<0)E20(0)|,0(t) 


-FM<0)E30(0>"'0(’> 


while 


(26) 


,2 

d m„ 


dx 


2 B10^E10^m0’ 


Since  E(0)B^q(0)  - 0,  2 (E(O)cIq)  = 0,  so  the  only  solution  which 

decays  to  zero  together  with  its  first  derivative  as  x ■+  °°  must  satisfy 


(27) 


E(0)mQ(x)  = 0. 


Furthermore,  multiplying  (26)  by  E^q(0),  provides  the  decaying  solution 


(28) 


-^E  (OB  (°)x 

E10(0)m0(T)  = e E10(0)m0(0), 


and  the  definition  (12)  of  E implies  that 


(29) 


-/E  (°)B  (0)x 

mQ(T)  = Q(0)e  XU  E10(0)m0(0), 


Finally,  the  differential  equation  for  p1()  has  the  decaying  solution 
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(30)  P10(t)  = - (/E10(0)B10(0))_1E10(0)ni0(T). 

Since  P^qCO)  and  E^p(0)m0(0)  are  arbitrary,  we  are  able  to  provide  m^ 
linearly  independent  solutions  to  (24)  by  specifying  either.  Higher  order 
coefficients  in  (22)  satisfy  nonhomogeneous  forms  of  (24)  with  successively 
known,  exponentially  decaying  terms.  The  decaying  solutions  (22)  are 
thereby  completely  determined  up  to  specification  of  either 


(31)  P1(0,u)  or  E1Q(0)m(0, u) . 


This  possibility  of  a choice  of  boundary  values  makes  the  boundary  layer  cor- 
rections with  0(p)  boundary  layers  more  flexible  than  those  with  0(e) 
layers  (cf.  (40)). 

The  classical  boundary  layer  correction 

oo 

(32)  (u2j,vA  l ) ^ l (u2j  ,y2£  ,p2t  ,P  )Mm 

1 L o ~ m lm  zm  ,5m 

m=0 


must  satisfy  the  system 


(33)  i 


= u2Aj  + y^l,  + y2B2i2  + B^., 

d^l  2 2 

— = Eri  + y Pjlj  + y F2^2  + F3S 

d£2  2 

~ e2j  + p F4«,1  + I-5«.2  + F6t3 


V di 


dK  P E3j  + p F7i1  + y Fg«,2  + F^. 


The  resulting  limiting  problem 
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dJf 


df'10 

= B,rt(0)a,n,  = E.n(0)jft  + F_n(0)l, 


dK  30v  ' 30’  dK  10  J0  30  30 


(34)  < 


di 


20 


da 


30 

dK  ~ E20(0)j0  + F50(0)5l20  + F60(0)2,30’  dK  = F90(0)''30 


has  the  unique  decaying  solution 


FqO(0>< 

So(K)  = e So(0)> 

jQ(K)  = B30(o)F90(0)a30(K) , 

(35)  1 a10(K)  = (e10(0)b30(0)f9q(0)  + F30(0)F9J(0))a30(K), 


and 


£20(k:) 


0CO  p (0)  (k-s) 

e 50  (E2Q(0)B30(0)F9J(0)  + F60(0))a30(s)ds 

K 


up  to  selection  of  £.^(0).  Higher  order  terms  satisfy  nonhomogeneous  forms 
of  (34)  and  are  completely  determined  up  to  the  initial  value 


(36)  a3(0,y). 

This  procedure,  then,  allows  m3  linearly  independent  boundary  layer  solu- 
tions (32)  to  be  constructed. 

The  terminal  boundary  layer  solutions  are  found  in  a manner  analogous 
to  those  at  t = 0.  Thus,  the  leading  terms  of 

U 

(37)  (pn,qj,pq2,uq3)  * Z 

m=0 


are 


“20(o)  ' -F50(l>  W°>- 


q30^°^  F90^^E30^^n0^0^  * 


"/E10(1)B10(1)° 

qio(a)  = e q10(0) 


nn(a)  = ®10<1)  ( l^i0("l)B10 (1)  n(°)  ' 


We  note  that  (12)  implies  that 


E(l)no(0)  = 0 


so  that  n^o)  = Q(l)E^Q(l)nQ(a)  and  these  leading  terms  are  therefore 
determined  up  to  the  initial  value  qin(0)  = (v/E10'(l)B  "0(1))"1E1  n(l)nn(0) 


'10'  '0 


Analogous  work  successively  provides  the  higher  order  terms  in  (37)  up  to 
specification  of  either 


q1(0,p)  or  E1Q(l)n(0,u) . 


Finally,  the  leading  terms  of  the  boundary  layer  correction 


(41)  (p2r,p2srs2,p2s3)  * £ ^V^lm’^m^Sm^ 

m=U 


are  given  by 


"F50(1)p 

S20(P)  = G S20(0) 


r0(P>  B20^F50^S20^ 
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sl0(p)  (E10(l)H20(l)F50(l)  + F20(l)F50(l))s20(p) 

f*  Fqq(1) (t-p)  _ 

s30(p)  - e (h30(l)B20(l)F50(l)  + Fg0(l))s20(t)dt 


and  the  solution  is  completely  specified  up  to  the  initial  m2  vector 


s2(0,u) . 


We  could  prove  the  asymptotic  validity  of  the  solutions  (20)  which 
we've  constructed  by  integral  equations  methods  (cf.  Harris  (1973)  or 
Vasil'eva  and  Butuzov  (1973)).  Although  that  proof  would  differ  somewhat 
from  the  classical  (regular)  ones,  we  regard  the  construction  of  the  solu- 
tions as  the  most  challenging  aspect  of  this  study  and  shall  not  further 
discuss  the  details  of  proof. 


4.  Fitting  the  Boundary  Conditions. 


Since  our  construction  of  the  outer  solution  (17)  and  of  the  thicker 
boundary  layer  corrections  (22)  and  (37)  distinguish  between  components  of 
Ex  and  E^x,  natura-*-  to  write 


x = + Qx2 


for  the  n and  m^  dimensional  vectors 


x3  = E(t)x(t,e)  and  x2  = E1Q(t)x(t,e) 


(We  experienced  a similar  separation  of  components  in  solving  the  reduced 
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9- 


and  u is  the  n + + m dimensional  vector 


(50) 


u(t)  = (x’  x'2  y | y'2  yp'. 


Furthermore,  the  expansions  (42)  imply  that  all  solutions  (20)  of  the  trans- 
formed problem  (6)  are  given  by 


(51) 


u ( t ) = U(t,M)k(p) 


where  the  square  matrix  U is 


1 

j x1(t,p) 

pm1(T,u) 

U2n1(a,y) 

h2r1(p,y) 

vj2i1('<>b) 

x2(t,p) 

m2(x,y) 

yn2(a,p) 

y2r2(p,y) 

V2J2<k  ,p) 

(52) 

U(t,vi)  = 

Y1(t,y) 

u pi(T’y) 

q1(a.w) 

y2s  (p,p) 

y2«.1(ic,y) 

1 

Y2(t,p) 

P2(t,m) 

pq9(cr,y) 

s2(p,p) 

2 

P &2 (k  > P) 

W3(t,p) 

P3(t,p) 

yq3(o,p) 

s3(p,p) 

^3(k,p) 

with 

fixed  boundary  values 

X1(j,y) 

= V j 

= 0 or  1, 

m2(0,u)  = 

Im1  = V°’v0 

(53) 

1 

1 

s2 (0, p) 

~ Im2* 

and  £3(0,p) 

1 


and  the  n + m + m vector  k(p)  is  partitioned  as 
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k(y)  = (k’(y)  y“k^(y)  y6k^(p)  yYk^(y)  y6k^(y))’ 


We  note  that  the  boundary  values  (53)  imply  that  p^(0)  = ~ (^E^gCOB^JCo) )~  , 
and  n2g(0)  = ( /E1q (1) B1Q (1) ) \ so  both  p^O.y)  and  n2(0,y)  are  non- 
singular for  e sufficiently  small.  Also  note  that  U(t,y)  would  be  a 
fundamental  matrix  if  the  corresponding  n + m + m^  dimensional  system  were 
homogeneous.  We  shall  select  a,B,Y,  and  6 so  that  the  k± (0) ’ s are 
0(1)  as  e ->  0 and  nonzero  if  the  particular  k^  is  not  identically  zero. 

Putting  (48)  and  (51)  together,  the  boundary  conditions  (2)  imply  the 
n + m linear  equations 


A1(y)k(y)  = c(e) 


for  the  unknowns  k^(y)  where 


(56)  A (u)  - Z (R . (y  ) S . (y2) )T( j )U( j , y) . 

j=0  J J 

An  additional  m^  boundary  conditions  result  if  we  impose  the  endcondition 

(57)  A2 (y)k(y)  = Elf)(j)x1(j  ,y2)  =0,  j = 0 or  1, 

required  by  consistency  with  (47).  (Note  that  has  rank  n^.)  Thus, 

the  boundary  value  problem  (1)  - (2)  will  have  a unique  asymptotic  solution 


of  the  form 


x(t,e)  \ / c (y 

= T(t)U(t,y)A_i(y) 
y(t,c)  / i o 


c(y2) 
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provided 


(H4)  the  (n  + m + m^)  x (n  + m + m^)  matrix  A(y)  = (A|(y)  A'(y))'  is 


nonsingular  for  e sufficiently  small. 


Because  the  boundary  layer  correction  terms  are  asymptotically  negli- 
gible away  from  one  endpoint,  we  can  considerably  simplify  calculation  of 
A(u).  Thus,  up  to  asymptotically  exponentially  small  terms. 


T(0)U(0, y)  = 


X1(0,y)  + Q(0)X2(0,u)  ym1(0,y)  + Q(0) 


OOy  j1(0, y)  + Q(0)j2(0,p)  \ 


3 3 2 

l P (0)Y  (0,y)  ~ P (0)Pl (0,y)  + Z P(0)p(0,y)  0 0 Z y2P,  (0)A.  (0,y)  + P,(0) 

£=1  M 1 £=2  * * k=l  K k J 


and  likewise  for  T(l)U(l,y).  These  imply  that 


(59)  A (y)  = (Ak£),  k = 1,  2;  A = 1,  ....  5 


where 


^ (R.(X  (j,y)  + Q(j)X„(j,y))  + S Z P.  (j)Y„ (j ,p) }, 

j=0  3 3 SL=1  11  * 


&12  * VPml(0,Vi)  + Q(0) ) + - S0(P1(0)p1(0,y)  + y E P^OJp^O.y)), 


A13  ^ R1^P  + PQ(1)n2(0,y))  + S^P^Dq^O.y)  + y Z Pfc(l)qj|  (0,y))  , 
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A14  ^ ^2R1(r1(0,u)  + Q(l)r2(0,u))  + S^p^Ds^O.p)  + 2 P^Ds^  (0,p)) , 

2 

A15  ^ Vu2jl(0,M)  + Q(0)j2(0,p))  + S0(p2  2 Pk(0).k(0,p)  + P (0)). 

k=l 


A2k  n,  E-^gCj  )X^(j  ,p)  with  j determined  in  (57) 


and 


(A22 ’ A23’  A24’  A25)  ^ { 


PElo(O)  (m1(0,p),  0,  0,  pj  <0,p))  if  j = 0 
and 

y2E10(D  (0,  n^(0,p) , r ^ (0 , p) , 0)  if  j = 1. 


Because  A^2  = 0(l/p)  (=  0(1)  only  if  S^P^(0)p^q(0)  = 0),  the  matrix  A(p) 
will  have  an  asymptotic  series  expansion 

oo 

(60)  A (p)  -V  - 2 S.pj. 

M j=0  J 

It  will  therefore  be  nonsingular  if  the  limiting  matrix  6^  is  nonsingular, 

although  A(p)  can  still  be  nonsingular  for  e f 0 if  6^  = 0.  If  6^  is 

-1  l-£ 

the  first  nonsingular  coefficient  in  (60),  A (p)  will  be  0(p  ),  so  the 

solution  (58)  of  the  given  problem  will  generally  be  unbounded  like  0(p  ). 

In  particular,  note  that  a bounded  solution  will  result  if  £ = 0 and  that 
the  powers  a,8,y,  and  6 in  (54)  are  integers.  Further,  the  limiting 
solution  within  (0,1)  depends  only  on  k (p)  , so  we  might  say  which  bound- 
ary conditions  are  appropriate  for  the  reduced  problem  (7)  (cf.  Harris  (1973))- 
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We  might  also  consider  the  possibility  of  nonunique  solutions  under  appro- 
priate orthogonality  assumptions  if  A(y)  is  singular. 

To  summarize  our  principal  results,  we  have 

Theorem: 

Under  hypotheses  (HI)  - (H4) , we  obtain  a unique  solution  (58)  of  the 
boundary  value  problem  (1)  - (2). 

5 . Natural  Boundary  Value  Problems, 
a)  Sample  Problem  1. 

Suppose  we  are  given  a problem  in  the  transformed  form  (6)  with  pre- 
scribed boundary  values 

(61)  x(0),  y^l),  y2(i),  and  /s  (!)  • 

Instead  of  actually  obtaining  A d(u),  we  can  apply  the  boundary  conditions 
in  (51)  to  obtain  a solution  with 


(62) 


k (0)  = E(0)x(0)  for  X.(0,p)  = I , 
i in 

a = 0,  k2(0)  = E1q(0)x(0)  - X2Q(0)k1(0), 

8=0,  k3(0)  = yx(l)  - Y10(l)k1(0), 

Y = 0,  k4(0)  = y2(l)  - ¥^(1)^(0), 

6 = 0,  k5(0)  = y3(0)  - ¥^(0)^(0)  - p30 


(0)k2(0). 


In  particular,  the  limiting  solution  within  (0,1)  will  satisfy  the  reduced 
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1 


problem  (7)  and  the  initial  condition  (16)  with  j = 0. 
b)  Sample  Problem  2. 

Suppose  our  problem  is  of  the  transformed  form  (6)  with  prescribed 
boundary  values 


(63)  x(l),  y1(0),  y2(l),  and  y3(0). 


Again,  the  unique  solution  is  readily  found  to  be  of  the  form  (51)  with 


(64)  / 


k1(0)  = E(l)x(l)  for  X1(l,p)  = In, 
a = 1,  k2(0)  = - /E10(0)Bl0(0)  (y1(0)  - Y10(0)k1(0)), 

B = "I.  k3(0)  = /E10(1)B10(1)  (E10(1)x(1)  - X20(l)k1(0)), 

Y = 0,  k4(0)  = y2(l)  - Y2Q(l)k1(0)  - g20(0)k3(0), 

6 = 0,  k5(0)  = y3(0)  - Y3Q(0)k1(0) 


and  the  limiting  solution  within  (0,1)  satisfies  the  system  (7)  and  the 
terminal  condition  (16)  with  j = 1. 
c)  The  General  Problem. 

As  our  special  problems  suggest,  we  can  seek  a solution  (51)  of  the 
transformed  system  (6)  plus  boundary  conditions.  That  problem  will  be 
uniquely  solvable  provided  the  reduced  system  (7)  can  be  uniquely  solved 
subject  to  appropriate  boundary  conditions.  We  must  suppose  that 
E1Q(0)x(0,e)  or  y^O.e),  E1q(1)x(1,e)  or  y^l.e),  y2(l,e),  and 

y3(l,e)  can  be  obtained  in  order  to  uniquely  determine  the  initial  values 
(31),  (36),  (40),  and  (43)  for  the  various  boundary  layer  corrections  of  the 
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solution  (20)  - (46).  In  order  to  solve  the  reduced  problem  (7),  separated 
boundary  values  for  E(j)X(j,e)  need  not  be  given  (as  in  the  sample  prob- 
lems). However,  if  the  prescribed  boundary  values  for  Ex  are  coupled  at 
t = 0 and  t = 1,  existence  of  the  solution  to  the  resulting  reduced  two 
point  problem  is  not  a priori  guaranteed.  One  must  always  be  able  to  solve 
the  given  conditions  for  y2(l,e)  and  y3(0,e).  Knowing  >’3(0,c),  for 
example,  one  would  solve 


y3 (0, c)  a,  Y3(0,p)k1(p)  + yap3(0,y)k2(p)  + y6k5(p) 


to  obtain  k^. 

Further,  it  is  essential  that  at  least  m^  + m3  boundary  values  be 
obtainable  at  t = 0 and  at  least  m^  + m2  boundary  values  be  obtainable 
at  t = 1 because  these  are  the  number  of  linearly  independent  boundary 
layer  corrections  decaying  at  those  endpoints.  In  particular,  we  cannot 
expect  to  asymptotically  solve  initial  value  problems  or  '-erminal  value 
problems  for  (1)  unless  we  artificially  restrict  boundary  values  to  appro- 
priate lower  dimensional  manifolds  (cf.  Hoppensteadt  (1971)). 

Since  the  solution  of  general  problems  (1)  - (2)  relates  crucially  to 
the  solution  of  simpler  transformed  problems  (like  our  sample  problems),  it 
is  convenient  to  solve  general  problems  in  terms  of  simpler  "natural"  ones. 
This  generalized  "shooting"  method  has  been  somewhat  developed  by  Keller 
and  White  (1975)  and  Ferguson  (1975). 
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